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Abstract 

We extend the Operator Product Expansion for Null Polygon Wilson loops to the Mason- 
Skinner- Caron-Huot super loop dual to non MHV gluon amplitudes. We explain how the 
known tree level amplitudes can be promoted into an infinite amount of data at any loop 
order in the OPE picture. As an application, we re-derive all one loop NMHV six gluon 
amplitudes by promoting their tree level expressions. We also present some new all loops 
predictions for these amplitudes. 
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1 Introduction 

Four dimensional gauge theories are, secretly, string theories [1]. This is really exciting both 
conceptually as well as pragmatically. In practice it means that if one is smart enough one 
can map four dimensional computations into two dimensional problems; these are of course 
way simpler. A remarkable property of some gauge theories is that they are integrable 
meaning that their two dimensional 't Hooft worldsheet is exactly solvable. The most known 
example is Af = 4 SYM also popularly known as the Harmonic Oscillator of Gauge Theories 
[2]. An object which has been under attack from several fronts is the planar S-matrix of this 
theory. One approach towards addressing the computation of scattering amplitudes using 
integrability, known as the Operator Product Expansion (OPE) [3], goes as follows 

• First we make use of a remarkable duality which tells us that planar scattering ampli- 
tudes are computed by the expectation value of (some supersymmetric) Null Polygon 
Wilson loops (NPWL) [4, 5, 6, 7]. Maximally Helicity Violating (MHV) amplitudes 
are computed by the usual bosonic Wilson loop [4]. Amplitudes with other helicities 
are given by a supersymmetric version of this loop [5, 6]. 

• Next, and most importantly, one needs to rephrase the problem of computing the four 
dimensional amplitude into a manifestly two dimensional language. To do so, we single 
out two of the edges of the polygon loop. We think of them as a quark/anti-quark pair 
moving at the speed of light and sourcing a color flux tube stretching between them. 
This color flux is thick in spacetime. However, it has an holographic description in 
terms of a two dimensional object [8]. That is the two dimensional description we were 
looking for. 
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Bottom 

Top 



Figure 1: The OPE picture. Flux tube excitations are created at the bottom and absorbed in the 
top. The flux is tick in spacetime but is 1+1 dimensional in AdS. Its excitations are integrable. 

• The choice of the two edges splits the polygon into a top and a bottom parts. We 
interpret the computation of the NPWL as describing the creation of excitations in the 
bottom part of the polygon which propagate in the flux tube and then absorbed by the 
top part of polygon, see figure 1. In other words we compute a transition amplitude in 
the two dimensional system. 

• The computation is now taking place in one (plus one) dimension. The energy of the 
flux tube excitations in planar Af = 4 SYM can be studied to all loop orders using 
integrability [9, 10, 11]. This means that the propagation part of the computation is 
under good control. The computation of the creation/absorption form factors within 
the integrability formalism is not yet known. A strategy to overcome this difficulty is 
to predict a big part of the result, coming from corrections to the propagation part in 
all possible channels. 1 That is, for all possible choices of pairs of edges. Then we try 
to bootstrap the full result by imposing crossing symmetry. 

This program was carried out in [14, 15, 16] for MHV amplitudes , that is for bosonic Wilson 
loops. Using this method, two loop results [17, 18, 19] obtained from direct computations 
were reproduced and infinitely many higher loop predictions were presented. One feature of 
all these OPE computations is that only single particle perturbations of the flux tube were 
considered. 2 To go beyond we need to better understand how multi-particle flux excitations 

^^More precisely, in the OPE approach, we compute the discontinuity of the Wilson loop in some conformal 
cross-ratios Uj. Similarly, the Unitarity method [12] for computing scattering amplitudes is based on the 
computation of discontinuities in Mandclstam invariants. It would be interesting to find a concrete relation 
between the two. This was partially achieved in the very interesting recent work [13] where discontinuities 
of Fcynman diagrams w.r.t. cross-ratios were taken. 

2 This resembles a lot the leading log expansion in Reggee theory and BFKL [20]. It would be very 
interesting to establish a more solid connection, see [21] for interesting advances in this direction. 
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One Loop MHV 



Tree Level NMHV 




Two Loop MHV One Loop NMHV Tree Level N 2 MHV 

Figure 2: One remarkable property of the Super Wilson Loop/Scattering Amplitude duality is that 
a N fc MHV scattering amplitude at I loops is given, roughly, by a k + / Wilson loop computation. 
This means that the known tree and one loop data from Scattering amplitude present an excellent 
window towards higher loop physics from the Wilson loop point of view. In particular, multi-particle 
flux tube excitations can in principle be studies with the already available data. 

behave. In principle, with integrability, all we need to know is how two particles interact. 
The scattering of N particles then factorizes in a sequence of (N — 1)! two body scattering 
events. From the MHV point of view, the critical missing point of data which, optimistically, 
should shed light over the full integrability structure is a three loop scattering amplitude. 3 
That seems like a very hardcore computation with the current Feynman diagramatic tools. 

In this paper we initiate the study of non MHV amplitudes from the OPE point of 
view. Our main motivation resides on a remarkable feature of the Wilson Loop/Scattering 
amplitude duality which is that a N fc MHV scattering amplitude at I loops is given by a 

3 The simplest of such amplitudes is probably the three loop eight point amplitude in a restricted kine- 
matics regime where all momenta of the external particles lie in the same plane [22]. 
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Wilson loop computation at, roughly, I + k loops [5, 6]. 45 A one loop MHV amplitude 
is given by a NPWL one loop computation with a gluon stretching from any edge to any 
other edge of the polygon. Similarly, some tree level NMHV component are given by the 
propagation of a scalar from one edge to another edge which clearly resembles the one loop 
computation. A N 2 MHV amplitude might be described by two such exchanges, resembling 
therefore a two loop computation and so on, see figure 2. Hence a plethora of multi-particle 
flux tube excitations abounds in the description of N fc MHV amplitudes, already at tree level 
and one loop! 

To be able to analyse this very exciting data through the OPE lens, we need to un- 
derstand what are the new ingredients in the OPE for supersymmetric NPWL. This is the 
purpose of the current paper. We will mostly focus on Next-to-Maximally Helicity Violating 
(NMHV) amplitudes; we will consider N fc MHV amplitudes and the multi-particle contribu- 
tions elsewhere. We will show how to bootstrap the full result for all six particle one loop 
NMHV amplitudes [25, 26, 27] from the tree level data [28]. We will also present some new 
all loop predictions for these amplitudes. 

This paper is organized as follows. In section 2 we will start by studying the tree level 
NMHV ratio function. We will explain how we can interpret known results from the OPE 
point of view. More precisely, we will re-write the results for six gluons as an expansion 
of particles flowing in the flux tube mentioned above. These expressions are the seeds of 
the higher loop computations. At one loop there are three types of corrections to the ratio 
function. Corrections to the energy of the particles, corrections to the form factor and a new 
two particle contribution, see figure 3. Due to integrability, the correction to the energy of 
the excitations is under full control. That correction captures what is known as the OPE 
discontinuity of the amplitude. We will be able to bootstrap the full one loop results from 
their discontinuities. To compute all discontinuities of the result, we will first compute a re- 
strictive set of discontinuities that have a clear OPE interpretation. All other discontinuities 
will be obtain from these restrictive set using SUSY Ward identities as explained in section 
3. In other words, the role of channel duality is played by the Ward identities. 6 In section 4 
we show how the OPE for NMHV loops can be used to bootstrap the one loop ratio function 

4 The fact that constant complexity slices are roughly the I + k = constant slices is also a feature of the 
all loop recursion relations [23]. In fact, it would be very interesting to understand how the OPE expansion 
hides inside the integrand. From the OPE point of view the Amplitude/ Wilson loop is divided into several 
different contributions with very precise physical meaning. It would be very interesting to try to isolate them 
at an earlier stage (already at the Integrand level if possible) hence simplifying their computation. 

5 The equivalence between the super Wilson loop and the scattering amplitude was proven at the level 
of the loop integrand [5, 6]. These observables however are UV/IR divergent and need regularization. A 
proper regularization of the super Wilson loop under which the duality with amplitudes stays untouched is 
not yet known; for example, dimensional regularization does not seem to work [7]. Here we assume that such 
a regularized object exists and we read the results from the known scattering amplitudes data. On the other 
hand, the null correlation function approach of [24] leads to basically the same super Wilson loop while their 
regularization is currently under better control. For us, this is not really important; all we care is that a 
picture with the two null lines sourcing the flux tube exists. Hence, the logic of the paper will not be effected 
is we replace any reference to the "Wilson loop" by "null correlation function". We thank E. Sokatchev for 
discussions on this point. 

6 One other nice outcome of the Ward identities is that they will allow us to bootstrap all one loop 
amplitudes. Curiously, this even works for components which vanish at tree level and hence start at one loop 
only. 
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Tree level Flux tube Two particle Form factor 

interaction contribution correction 



Figure 3: T]ii]^r\^r]% component at tree level (left) and at one loop level (right). At one loop there 
are three contributions that kick in as indicated in the figure. From the OPE point of view the 
interaction with the flux tube is under control since we can compute the energy of the flux tube 
excitations exactly using integrability [9]. 

and in section 5 we present some new all loop predictions for the NMHV ratio function. 



2 Tree Level 

The observable we will study in this paper is the ratio function for NMHV amplitudes [29]. 
In this section we study this observable at tree level which as we explain below is the seed 
to the higher loops. It is obtained by extracting components from the NMHV generating 
functional [29, 30, 31] 

A NMHV tree 

Kr HVa °°=%^= E [i,m+w+2] (i) 

n l<i<j<n-2 

where the five bracket is the so called R-invariant given by 

[abode] — ^ ^ a ^ cc ^ e ) + Vb{cdea) + rj c {deab) + r) d {eabc) + rj e (abcd}) 
' ' ' ' (abed) (bade) (cdea) (deab) (eabc) 

and (abed) = eabcdZ^Z^ is the determinant of four dual momentum twistors. We 

assume the reader is familiar with these formulae (for nice reviews see [32, 27, 33, 34] and 
references therein). The main purpose of writing them down was to set the notation used 
throughout the paper. Lorentz, spinor and R-charge indices are always omitted unless their 
placement is not obvious. As a final disclaimer, we also assume acquaintance with the basic 
OPE ideas spelled out in [3], see also [14, 15, 16]. 7 In this paper we will focus mostly on the 



7 In particular, see [3] for more details on the concepts of flux tube, flux tube excitations and their 
anomalous dimensions, reference square, R T x R a x SO(2),p symmetries of the reference square, OPE channel 
defined by two null lines, SL(2) T symmetries of the two null lines and also meaning of the parametrization 
variables T,a,4> and their conjugate charges E 7 p,m. For the idea of the SL{2) T Casimir operator acting 
on part of a null polygon see sections 3.5 and 5.2 in [15] and also [16]. For more details on the flux tube 
anomalous dimensions see [9, 14, 15]. 
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new features of the OPE which arise when considering supersymmetric loops. 

For simplicity, in this paper, we will consider NMHV amplitudes with n = 6 particles. 
That is we will consider the Super Null Hexagon Wilson Loop. The generalization to more 
particles can be done following [16]. 

The OPE description involves thinking of particles being produced in a region of the 
polygon, propagating in a gluonic flux tube and being absorbed at some other region of the 
polygon. Hence, most symmetries of the scattering ampltiudes, in particular ciclicity and 
supersymmetry, are not manifest. In fact, when analyzing amplitudes from the OPE point 
of view, we found it sensible to pick components instead of considering the whole generating 
function (1) at once. Let us start with a simple example, the rj2rj 3 r] 5 r]Q component 

^NMHVtreel _ ^(2356) _ \ /g\ 

6 \mmmv6 ~ ~ (2356) ' 

and analyze it from the OPE channel where the two null lines of the expansion coincide with 
edge 1 and 4 (see figure 3. a). 

An important new conceptual novelty one encounters when dealing with objects like (3) 
is that these carry helicity weights and are therefore not pure functions of conformal cross 
ratios. That is, the momentum twistors live in projective space and (3) is not a c-number. It 
only becomes a number, measurable in a collider, when contracted with the corresponding 
polarization vectors. 

To apply the OPE procedure to these helicity carrying objects, we will first fix an helicity 
gauge by picking some momentum twistors Z a for some initial reference hexagon. We then 
generate a three parameter family of momentum twistors Z(t, a, 0) by acting with the three 
conformal symmetries of the reference square R T x R a x SO (2)^ on the bottom twistors, that 
is 

M(r, cr, <p)-Z a , for a = 2,3 ^ 
Z a , for a 7^ 2, 3 

where M(r, a, <fi) is the R T x R a x SO{2)^ transformation mentioned above. 8 For the choice 
of reference square used in [15], see also appendix A, we find 

n^[Zi] = 4 [ cosh(a) C o S h(r) + cos(0)] (5) 

The overall normalization of the right hand side of (5) depends on our choice of helicity gauge 
Z a , see (23). We would like to express (5) as a sum of flux tube excitations propagating from 
the bottom to the top. These can be classified into SL(2) primaries and descendants. The 
NMHV amplitude (3) is obtained from a scalar propagator between cusp Z 2 AZ S and Z 5 AZ 6 . 
Therefore, all primaries are expected to be scalars. All such primaries are 

d™<j> AB , df(P AB , where m = 0, 1, 2, . . . (6) 

8 To be more precise, one acts on the bottom cusps. These are located at Z\ A A Z3 and Z3 A Z±. 

That is, they also depends on Z\ and Z^. However, Z\ and Z4 only transform under R T x R a x SO{2) t f ) by 
an overall rescaling (see (23)). For the channel defined by edges 1 and 4, we will only consider those NMHV 
amplitudes that do not carry 771 and r/4 helicity weights. These are invariant under overall rescalings of Z\ 
and Z4. 
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Here, d Z} dz are spacetime derivatives in the transverse directions to the reference square, 
carrying charge ±1 under SO {2)$. There is a very simple check one can do to confirm that 
indeed, (5) has a decomposition in these excitations and their conformal descendents. The 
twist of these excitations is E = A — 5 = 1 + \m\. The corresponding SL(2) Casimir is 
C = E(E — 2) = (m 2 — l). 9 In particular, it is independent of the sign of m and therefore its 
action on excitations (6) is given by the differential operator — (9| + 1). On the other hand, 
the SL{2) Casimir can be represented as a second order differential operator 

— J ^ ~\~ J 2 J ^ 

acting on the bottom of the polygon. Its action on the ratio function can be computed as 10 

3 

C o Ufa] = ^(-l) a+1 d 2 e U [Z lisU , {Z top }, Z MU {e i9J « o Z bot }} \ g=Q 
a=i 

We conclude that the ratio function has the expected decomposition if and only if it is 
annihilated by the second order differential operator 

P = C + (<9j + 1) . 

Indeed, we find 11 

P o n {2356) = . 

This is first non-trivial evidence for the OPE expansion of super loop! Let us check another 
component where we also expect scalar excitations to be flowing, 

^(2256) [Z j (3561) csch(r)e^ + coth(r) 



(2356) (1256) 4 cosh(o-) cosh(r) + 4 cosh(</>) 
Once again we find 

P o ft (2256) [Z,] = 0. 

It is useful to convert the action of C on TZ^^ [Zi(r, a, 0)] into a differential operator C^ J ' fe ^ 
acting on a, r, 0. We find 

C (2356) = [^ + 2c oth(2r)a r + sech 2 (r)^] (7) 
C (2256) = ^2 + 2coth(2r)9 r + sech 2 (r)<9 2 ] + [sech 2 (r)(l - 2d a ) - csch 2 (r)] (8) 

Note that the operator C, when converted into a differential operator acting on a, r, <p, takes 
different forms when acting on the ^2^3%% component (7) and on the ^2 ? ?2 ? ?5 ? 76 component 
(8). This arises because we are dealing with functions carrying helicity weight. The reader 
who is not puzzled with this statement should probably skip the next paragraph (in smaller 
font) in a first reading. 



9 Our convention for the SL(2) Casimir here differ by a factor of 4 with respect to [14, 15]. 

10 In the paremetrization presented in the appendix, J\ = I2X2 ®<Ji, J2 = i^2x2 ® &2, J3 = ^2x2 ® o"3 where 
the cr a 's are Pauli matrices. Note that ^ r i g ht and Z\ e ft are invariant under the SL(2) T action and therefore 
the result is the same whether we act on them or not. 

11 See also [34] 
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To understand why the differential operator C^ kl ^ takes different forms acting on different ratio 
function components lets us be more general and consider an arbitrary conformal invariant function 
[z a ] carrying helicity wight r\ir\j"(]k°f]i- Ratio function components are examples of such functions. 
Even though nS l i kl ) [Z a ] is conformal invariant, it is not a function of the conformal cross ratios only. 
Instead, it transforms under overall resealing of the twistors {Zi, Zj, Zf., Z{\. The ratio = 
%( l 3 kl )[Z]/lZ tyl i k ^[Z] however, is independent of the helicity gauge chosen and is therefore a function 
of the conformal cross ratios u, only. In other words, 

n^ kl) [z a ] = h ( - ijkl \u 1 ,u 2 ,u 3 )n^{z a ] (9) 

where (4) and tij's are the cross ratios computed from the transformed Z's (24). The parametrization 
of Mi, U2 and 113 in terms of r, a and <f> is given in the Appendix, see (25). (Implicitly, we assumed that 
the tree level ratio function component TZ^i kl ) ^ 0; if TZ^ kl ^ = one can use an arbitrary reference 
conformal invariant funcrion carrying (ijkl) weight.) Under an SL(2) transformation g acting on the 
bottom twistors, 

H dm [Z ] -> fe<««) (g or, go a, go <^« w > [Z 1)9 ■ Z 2 ,g ■ Z 3 ,Z 4 ,Z 5 , Z 6 ] (10) 

where jot and g o a stands for the action on r and a dictated from the transformation of the cross- 
ratios (26), see appendix. This action is of course universal. However, the last term in (10) has helicity 
weight and transforms differently for different components. That transformation gives a non trivial 
contribution to the differential operator CS % i kl ) and is the source for the difference between C*- 2356 ' 
and C^ 2256 ) for example. A more thorough discussion of this point is given in appendix A. 2. 



We will now decompose (5) in the scalar excitations propagating from the bottom to 
the top. The propagation of a conformal primary and its conformal descendants are conve- 
niently packed together in a conformal block. These are functions of the form (9) r Hy^ kl > = 



e %pa J-Ep(r) that solve the differential equation 



12 



(ijkl) 



a ipcr 



=E{E-2) 



E,p 



e ipajr 



(ijkl) 



E,p 



T 



The relevant solution of that equation with (ijkl) = (2356) is given in [15] and reads 



4TV)=sechV) .Fx 



^,^Usech 2 (r) 



We find that 

00 

^(2356) [z]= 

m=— 00 

where the form factors are given by 

4 2356) (P) = \(-V m B 



dp im<f>-ipa r (2356) T (2356) , 
271- m ^• r |m|+l,pV 



m 


+ 1 + ip 


m 


+ 1 - ip\ 




2 




2 J 



12 We use a slightly different notation here compared to [15, 16]. Namely, J 7 ^p l \ T ) 

The conformal cross ratio <fi do not enter the equation because it does not transform under the SL(2) 
transformation g o <p = cf>. 



11) 



12) 



(13) 



(14) 
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V2V2V5V5 



Figure 4: Some components with a good OPE expansion in the channel defined by the null edges 1 
and 4. The first two are given by scalars propagating from bottom to top and the third is described 
by propagation of fermions. 



For our other example we find exactly the same expansion with extremely similar blocks and 
form factors, (the two minor differences are colored) 



T^ p 56 \r) = tanh(r)sech s (r) 2 F X 

e 56) w = e 56) (P-2^ m >o). 



^ E + f + 2 ,E,sechV) 



It is very curious to notice that these expansions are shockingly similar to each other and 
also to the expansions arising for the MHV hexagon; compare for example (5), (13) with 
equations (75), (76) from [15]. The form factors (14) are exactly the same as those arising in 
the expansion of Uruny This seems to indicate that there is some sort of hidden symmetry 
relating them. In fact, this statement is even more reinforced when we consider other tree 
level components with a good OPE expansion in the channel defined by the edges 1 and 4, see 
figure 2 for examples. In terms of momentum twistors these component look very different 
from each other. However, when re-written as a decomposition over flux tube excitations, 
we found that they all look roughly the same up to very minor redefinitions of the conformal 
blocks and form factors which typically just involve some simple shifts in p and ml These 
similarities are partially a consequence of supersymmetry which tells us that only five NMHV 
components are actually independent. All other components can be obtained from these 
ones through SUSY Ward identities [35]. There ought to be a more supersymmetric way of 
packaging all these expansions together into some sort of SUSY blocks and/or SUSY form 
factors. 
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Figure 5: Ward Identities allow us to relate any six point scattering amplitude to a base of five 
components. In particular this allows us to access the behavior of any amplitude in any OPE 
channel. In the figure the component which clearly does not have a natural expansion 

in the channel defined by the null edges 1 and 4 is expanded in terms of five components that do 
admit a neat expansion in this channel. This can be done at any loop order. Ward identities replace 
channel duality in this problem. 

3 SUSY 

In this section we will review how different components of the NMHV ratio function are 
related to each other through supersymmetry [35]. We will also study some implications of 
cyclicity and parity. 

As explained in the previous section, when considering the OPE expansion in some chan- 
nel, there are particularly nice components which admit a simple expansion in that channel. 
For example, the amplitudes in figure 4 have a nice expansion in the OPE channel defined 
by the null edges 1 and 4. Naively, these components are tricky to expand in other chan- 
nels. On the other hand, to bootstrap the amplitudes at higher loops we need to be able to 
expand any component in any possible channel! SUSY Ward identities provide a neat way 
out. Using them we can rotate between components are infer the OPE expansion of some 
components in channels where naively that OPE expansion is hard to predict directly, see 
figure 5 for one such example. In other words, they replace/implement channel duality. This 
will be crucial when bootstraping one loop amplitudes from their OPE discontinuities in all 
possible channels. 
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The generating functional of NMHV ratio functions is superconformal and dual super- 
conformal. Hence it ought to be given by a linear combination of R-invariants (2). For the 
hexagon only five are independent [35]. Without loss of generality we can write 



n mwmv = [6 2 34 5] 21+ [612 35] 93 + [6124 5] £ + [61234]£ + [61345] £. (15) 

Next, we express 21, . . . , (E in terms of five components that we choose. For example, a 
particularly nice choice of components is the one in figure 5, 

7^(2356) .^(2256) ^(3356) ^(2355) ^(2366) Qg\ 

For example 

(1235) (23 45) (2355) (1236) (23 45) (2356) 
(1234) (1234) 

and so on (see appendix for explicit expressions). Note that these expressions are valid at 
any loop order, they are simply encoding SUSY Ward identities [35]. Actually, all we need 
to compute is two components!, for example the two components 7^. < - 2356 - ) ,7£( 2256 * ) considered 
in the previous section are enough. The other three of our base can then be obtained by 
simple symmetries of the hexagon, 

.^(2256) _ -^(3356)1 _ _ ^(2355)1 _ _ ^(2366)1 /yj\ 

I cr— > — (J ' I c/> — > — <6 I <fi— > — (ft , cr— > — a ' ^ ' 

We know how to promote all five amplitudes (16) in the OPE channel defined by the null 
edges 1 and 4. This was the reason why we chose this particular base. On the other hand, 
some other components such as 7?/ 1144 ) do not have an obvious expansion in this channel, 
see figure 5. However, given (15), we can now write any component as a linear combination 
of the base of amplitudes (16). In this way we can study the behavior of 7^ 1144 ) in a channel 
that was naively inaccessible, see figure 5. 

To summarize: we can study any channel we want, for any component we want, 

by choosing a base of five components with a good OPE expansion in that channel and by 
writing TZS 1 ^ in terms of that base. 

Using these tricks we can easily compute the one loop discontinuities of all NMHV com- 
ponents in all possible channels. With this information we can bootstrap the full one loop 
ratio function. This is the subject of the next section. 



4 One Loop 

One loop amplitudes will have discontinuities at u a = 0. In the Euclidean sheet we do not 
expect any other discontinuities [15]. The three discontinuities D a can be computed by the 
OPE in the three possible channels [15]. For MHV this suffices to determine the full two 
loop remainder function. In this section we show that the same is true for the NMHV ratio 
function at one loop. 
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From Tree Level to One Loop OPE Discontinuities 



Consider for example the two components studied in section 2. Their OPE tree level decom- 
position is given in formulae (13) and (15). We obtain the discontinuity in the u 2 = cosh" r 
channel by simply dressing these expansions by the anomalous dimensions of the scalar 
primaries flowing in the tube [15]. That is 

D (i m = J2 J ^ Cgf"> (p) ^2,+Ur) 7i +H (p) , 



where (ijkl) € {(2356), (2256)} and [9, 15, 14, 10] 

1 + \m\ + ip 



7i+|m](p) = 2# 



1 + \m\ — ip 



2^(1) 



Preforming these sums we get the corresponding OPE discontinuities of the one loop results 



D 



(2356) 



K 



D 



(2256) 



(2356) 2 , 
tree ^ lu 6 

log 



^,(2256) „ 
'nree z 



l-u 2 
"1 



U 2 



U 3 log(li 3 ) "I - Ml + M 2 + M 3 + y/UiU 2 U 3 \jl + ± 
1-W3 +1 - Ml + M2 + M3 + V /M 1 M 2 M 3 (" ~ ^ 



Here, (« = e J<?i ) and we have factored out the tree level results in order to express the result as 
a function of conformal cross ratios only, without the helicity weights ambiguities discussed 
in section 2. We can now get the other three components in (16) from simple parity and 
ciclicity symmetries (17) 



D 



(3356) 



D 



(2355) 



D 



(2366) 



^(2256) 9 
'^tree z 



-r>(2256) 
'Mree z 



-r>(2256) 9 



log 



log 



log 



U 3 



1 - U 2 



log(tti) -1 - "3 +«2 + «1 + ^U!U 2 U 3 Ut + J 



«1 



1 - u 2 



U3_ 

- u 2 



1 «1 +1 - U 3 + « 2 + «i + y/UlU 2 U 3 (u - i 
M 3 log(u 3 ) -1 - «1 +«2 + «3 - v/UlM 2 U 3 (/i + J 



1 U 3 +1 - + M 2 + M 3 + y/UiU 2 U 3 \J1 - ij 

Ml lQg(m) -1 ~ % + M2 + Ml - y^l^3 (/i + 
1-Wl +1 - M3 + M 2 + Mi + y/U^U^ \ji - 



Discontinuities from Ward identities, Ciclicity and Parity 

We have just computed discontinuities of the chosen base of five amplitudes (16) around 
u 2 = 0. Using the decomposition of the ratio function in terms of these five amplitudes (35) 
we can now read the discontinuity of any component in this channel. We simply need to put 
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together the five discontinuities computed above with the correct prefactors as predicted by 
(15). For example, for the example in figure 5 we find 



D 



(1144) 



K 



(1144) 
tree 



-21og(l-u 2 ) + 



2u 2 (ui - u 3 ) \og(u 3 /ui] 
(1 - u 2 )(l -u 1 - u 3 ) 



This was one of the tricky components from the u 2 = OPE channel point of view, see 
figure 5. Using cyclicity of the generating function, we can now read the discontinuities of 
any amplitude in the other two channels as well. These are the discontinuities around U\ — 
and m 3 = 0. For example, 



(2356) 



D 



(1245) 



7^(2356) ^(1245) J 



D 



(2356) 



D 



(1245) 



7^(2356) ^(1245) J 



Ui—>U i + 1 
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Bootstrapping the Full Amplitude From its Discontinuities 



It is now a very simple task to reconstruct the full function for our base from those disconti- 
nuities. One scientific way to proceed would be to use symbols [19, 15]. This works perfectly. 
On the other hand, we are dealing with a very simple expressions of transcendentally degree 
two. Hence, will proceed in a more pedestrian and less fancy way and illustrate, on a simple 
example, how to obtain the form of the full result from the knowledge of its three OPE 
discontinuities. 

Let us consider the particularly simple component 7?/ 2356 ); the other cases are roughly 
of the same complexity. Its three discontinuities can be derived from the OPE as explained 
above. One obtains 



(2356) 



D 



D 



(2356) 



(2356) 



2K 



2K 



2TZ 



(2356) 
tree 

(2356) 
tree 

(2356) 
tree 



logM 2 - logU 3 - log(l - Ut) 
log Mi + logM 3 - log(l - U 2 ) 
\ogU 2 - log Mi - log(l - U 3 ) 



(19) 
(20) 
(21) 



Consider first the discontinuity (19) in the channel u\ = 0. What function has logu 2 
logM 3 — log(l — Mi) as a discontinuity? There are two obvious options: 



log Mi ( log M 2 - log M 3 - log( 1 - Mi 



+. . . 



or 



logMi(logw 2 -logM 3 ) -Li 2 (l-Mi) + . 



13 Note that in some cases like the component (3556), the OPE discontinuity is associated with the corrected 
propagation of a fermion instead of a scalar, see figure 4. The dependence of the fermion anomalous dimension 
on the SO{2) ( f, charge m is different from the scalar one by a 1/2 shift: 7i+| m +i/2| instead of 7i+| m |. Still, 
the discontinuities for the two cases are simply related to each other. The way it comes about is that in the 
fermionic case, the charge m runs over half integer values so after factoring out an e - ^/ 2 , the dependence 
on the SO{2) ( j > charge is the same. 
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where . . . stands for extra terms without any discontinuity around U\ = 0. It is clear that 
we should pick the second option because the first option has a discontinuity at U\ — 1 and 
the only allowed discontinuities in the Euclidean sheet are at u a = and u a = oo, when two 
cusps become null separated. Hence 



,£,(2356) o (2356) 

'^onc loop — ^ 'Mree 



log Mi \ogu 2 - log Mi iogu 3 - Li 2 (l - Mi) + . . . 



(22) 



Next we look at the second discontinuity, in the channel m 2 = 0. The first term in (20) is 
already accounted by the first term in (22); the second term leads us an extra logM 2 logM 3 ; 
the last term gives rise to another Li 2 by the same reason as argued above. Hence 



--,(2356) r>-i-,(2356) 
'"^one loop — ^ 'nree 



logM 2 log(MiM 3 ) - loguilogM 3 - Li 2 (l - Ui) - Li 2 (l - M 2 ) + 



We only need to consider the last discontinuity (21). We easily see that to reproduce it we 
only need to add another polylogarithm to the result. That is, 



--,(2356) <-,--, (2356) 

"-"one loop — ^ 'Mree 



logU 2 log(MiM 3 ) - logUilogM 3 - ^Li 2 (l - U a ) + 



a=l 



where c is a constant of transcendentally degree two which we will soon fix using collinear 
limits. This is the prediction of the OPE approach. 11 We can now compare it with the 
known results computed in the literature by more direct methods [25, 26, 27]. We find a 
perfect match with the OPE prediction. 

All other components can be derived in the same way. Due to the Ward identities reviewed 
in the previous section, all we need is four more components: with a base of five amplitudes 
we can extract any component through (15). Furthermore, as also explained above, only 
one out of these four is really new, the other ones are trivially obtained from this one (17). 
Proceeding as above we derive that missing component, 



,£,(2256) _ r> <i-,(2256) 
'^one loop — * '^-tree 



log(Mi)log(M 2 ) + C ' + ^Li 2 (l 



a=l 



+ 



y>2 



ui+u 2 +u 3 - 2miu 3 - 1 + (/! + -) Juiu 2 uA log(u 3 ) log 

2(1-Mi)(l-M 3 ) 



We are done concerning one-loop except for the two constants c and d which we did not fix 
yet. These can be easily fixed from collinear limits to be c = d = — 7r 2 /3. 



The idea is that many components like ^on^ioop/^tree^ or ^on^itop/^tree 34 '' ou ght to vanish 
in the r — > oo limit. That is, any component that do not involve the edges {2,3} or the 
edges {5,6} will approach the NMHV pentagon in the r — > oo collinear limit. The NMHV 
pentagon is an MVH amplitude so that the ratio function does not receive loop corrections. 

14 Implicitly we also assume that at a given loop order the transcendentality degree of the result is uniform. 
Empirically, this seems to be a true property of AT — 4 SYM. 
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These components can be written in terms of the amplitudes derived above. In this way one 
concludes that 




This concludes the derivation of all one loop NMHV six point amplitudes from the OPE 
point of view. Note in particular that using the Ward identities we also get those one loop 
NMHV amplitudes that are zero at three level. Next we move to higher loops. 



5 All Loop Predictions 

Using the OPE expansion we can easily predict an infinite amount of data for higher loop 
NMHV amplitudes. For example, at n loops we expect [3] 

«S = A (n,C Vi, «2, «a) + A (n,1 W«2, us) log(«0 + ■ • • + ' n) K «2, «a) log(«i) B 

where D^' k \u a ) have regular expansions in powers of U{ around Ui = 0. Normally each 
jj{n,k) sYiovld have an extra index (2356) which we omitted to render the expression lighter. 
In the previous section we showed how to compute the discontinuities Dj- 1 ' = D% in all 
possible channels and we saw that they are sufficient to reconstruct the full NMHV result at 
one loop. The highest discontinuity D^ n ' n ^ at n loops is equally simple to get. For example 
[3] 

m=—oo 

where J 7 , C and 7 are given in (12), (14) and (18). To compute the maximal discontinuity in 
other channels we use the Ward identities exactly as above. The maximal discontinuities of 
other components are computed in the same way. 

This contribution is the easiest to extract from the OPE point of view. To bootstrap 
the full two loop NMHV ratio function one also needs to compute D\ . That requires a 
better handling of multi-particle flux tube excitations and is therefore beyond the scope of 
the present paper. 
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A Details on the Hexagon Computation 

A.l Parametrization and cross-ratios 

As mentioned in the text a family of hexagons is generated by starting with some random 
twistors Zi and acting with the symmetries of a reference square on a subset of those twistors, 
those constituting the bottom part of the polygon. A nice parametrization which we use 
throughout the paper is 



(z 1 \ 


( 


zA 


z 2 


M 


z 2 


z 3 


M 


Zz 


z 4 




z± 


z 5 




Zs 




\ 


zj 



( 



1 
1 

V o 






1 





-e 2 
-e 2 



1 

1 



. \ 





1 



(23) 



■1 / 



The left and right null lines correspond to {Z\ cit = Z 1: Z\ cit = Z§ A Zx A Z 2 } and {Z Tight = 
Zi, bright = Z 3 A Z4 A Z 5 }. The top and bottom twistors of the reference square are Z top = 
(0,0,1,0) and ^bottom = (0,0,0,1). The SL(2) group which preserves the two null lines 
{Zieft, Zi e ft} and {^ r i g ht) bright} acts on the two last slots of the twistors as I 2X 2 ® 9sl(2)i see 
footnote 10. A base of cross ratios 



"1 



u 2 



reads [15] 



(3456) (6123) 
(3461) (5623) ' 

e a sinh(r) tanh(r 



(1234) (4561) 
(1245) (3461) 



U3 



u 2 



2 cosh(cr) cosh(r) + 2 cosh(0) 
In particular, as an expansion at large r, we have 



cosh (r) 



log Mi = y^ci, w e~ 



logw 2 = -2r + c 2 , n e" 



n=0 



n=0 



(2345) (5612) 
(2356) (4512) 



u 3 = 



logw 3 = ^c 3jn e" 

n=0 



(24) 



(25) 



so that computing the term linear in r is the same (up to a factor of —2) as computing the 
discontinuity in u 2 [15]. Finally, it is also useful to have an inverse map relating directly r, 
a and 6 to cross-ratios, 



■2<r 



(1236) (1245) (3456) 
(1256) (1346) (2345) 



sinh (r) 



(3451) (4612) 
(3412) (4561) 



(1236) (1456) (2345) 
(1234) (1256) (3456) 



(26) 
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A. 2 Details on the Casimir operator(s) 

In this section we give a detailed derivation of the SL{2) Casimir differential operator C^ w ) 
acting on a conformal invariant function "H^ 7 ^ carrying helicity weight 77i%%7# (9). For 
concreteness, we will follow the example where (ijkl) = (2256) and C*- 2256 ^ takes the form 
(8). The generalization ot other components will be clear from the text. 

The tree level ratio function for the component (2256) is 

^ (2256) lz i = (3561) = _ csch(r)e^ + coth(r) 

[ ll (2356) (1256) 4 cosh (a) cosh(r) + 4 cosh(0) 1 ' 

The SL(2) transformation which leaves the null lines associated to the twistors Z\ and Z^ 
invariant reads 

/ 1 \ 
10 
0a b 
\ c / 

We then act with this transformation on the bottom twistors Z2 and Z%. We use the notation 



9 



(28) 



(1256) 9 = det{Zt,gZ 2 , Z 5 , Z 6 } (29) 
etc. Then we define the action of SL(2) on r, o and (f) using (26). For example 



, (1236) 8 (1245) a (3456) ! , 1 + (a - c)(ae^ - 6) 

g ° " = '° S » / <1256),<1346).<2345), = " + '° S V 1 + (« + c)(«» + 6) (30) 

and so on. In particularly, we find that g o cj> = 0. Two important quantities are 

- ( J 6 tTl256) g <*> 



csch(g o r)e 9O<T+0 + coth(p o r) 
4 cosh(g o a) cosh(g o r) + 4 cosh(0) 



Note that, despite (27), these two quantities are not the same! This is the most important 
point. Their ratio is a nontrivial function of a, b, c. Then we consider a general conformal 
invariant function "H( 2256 ) of the Zj's carrying helicity weight (2256). For the family of 
twistors {Z 2t 3 = M(r, a, 4>) ■ Z 2j3 , Z a = Z a for a ^ 2,3} it can be represented by a function 
of the three parameters r, a and 0, %( 2256 )(r, er, 0). Under the SL(2) transformation (29) of 
Z 2 and Z 3 , it transforms as 

-1/(2256)/ /\ 0/(2256)/ j. \ 7o 2256) [Zi, pZ 2 , #Z 3 , ^4, Z 5 , Z 6 ] 
0O?^ j T,(T,0 )=?^ T ,9 W^?2256W? TyH 32 

The second factor on the right hand side captures the transformation of the helicity factor. 
Now, let gg be the SL(2) transformation (29) with 



a b 

c ^tl 



(33) 



where <7j are the Pauli matrices. The action of the SL(2) Casimir on "H( 2256 ) reads 



C (2256) Q ^(2256) = ^ Q 2 g& Q ^(2256) 



a=l 



(34) 



(9=0 



([<9 2 + 2coth(2r)9 r + sech 2 (r)<9 2 ] + [sech 2 (r)(l - 2d„) - csch 2 (r)]) H (2256) 



For other components, the action of the Casimir is obtained is the same way. Since the 
transformation of TZy^ kl i is different for different components, so is the differential operator 
<C^ kl \ For the 2356 components for example, it turns out that the total helicity weight 
is invariant under the SL(2) action. As a result, if( 2356 ) transforms as a function of the 
conformal cross ratios. This is why the differential operator C^ 2356 - 1 (7) is exactly the same 
as the one we found for the hexagon reminder function in [15]. 



A. 3 Explicit solution to the Ward identities 



Ward identities, studied in [35] and in section 3, allow us to write the NMHV functional in 
terms of five components alone. This is extremely useful for the OPE program. It allow us to 
study any component in any channel by choosing appropriately the base of five components. 
For example, for the choice (16) we have 



NMHV 



= [6,1,2,3,5] 
- [6,2,3,4,5] 
+ [6,1,2,4,5] 



TZ {2355) (1235) (2345) + ft (2356) (1235) (23 46)" 




(1234) 

( 



' (2356) / (2346) (2351) (5612) (5643) \ 
(1234) (1456) J 



n (2256) ( 1256 ) (2456) 7^ 2355 > (1235) (2345) 



(1456) 



(1234) 



- [6,1,2,3,4] 



(2356) / (1236) (2345) (1235) (2346) \ 
1 (1234) (1234) 1 



) 



^ft( 2 355) (1235) (2345) ft( 2366 )( 1236) (2346) 



(1234) 



(1234) 



+ [6,1,3,4,5] 



(2356) / (5612) (5643) _ (2345) (2361) \ 
1 (1456) (1234) ' 



) 



(1235) (2345) 7^ 3356 ) (1356) (3456) 



(1234) 



(1456) 



(35) 
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